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INTRODUCTION AND SlfrMARY 

This report is a summary of research activities performed under grant 
NAG-I-6. The objective of the research has been to study the design of 
adapt ive/leaming control systems for the control of large flexible structures. 

In 1 an adaptive/leaming control methodology for flexible space structures 
was described. The approach was based on using a modal model of the flexible 
structure dynamics and an output-error identification scheme to identify 
modal parameters. The identification scheme was tested on simulated data 
characteristic of the solar electric propulsion (SEP) array. In 2 a least- 
squares identification scheme was proposed for estimating both modal parameters 
and mod al-to- actuator and modal -to-sensor shape functions. The technique was 
applied to experimental, data obtained from the NASA Langley beam experiment. 

In 3 a separable nonlinear least-squares approach was described for estimating 
the number of excited modes, shape functions, modal parameters, and modal 
amplitude and velocity time functions for a flexible structure. A digital 
computer program was delivered to the NASA- Lang ley Research Center and was used to 
process experimental flexible beam data to obtain estimates of the parameters 
and functions mentioned above. Because least-squares residuals are computed 
during the identification procedure , the method provides its own monitoring 
of identification quality. In 4 a dual -adaptive control strategy was developed 
for regulating the modal dynamics and identifying modal parameters of a flexible 
structure. It is well known -that input signals that are optimum for identifying 
system parameters often yield a dynamic response characterized by large ampli- 
tude modal variations. Hence a min -max approach was used for finding an input 
to provide modal parameter identification while not exceeding reasonable bounds 
on modal displacement. The approach was tested using simulated beam data. 


FUTURE PLANS 


Our future research plans include the following: (1) simplification of 
the separable nonlinear least squares corputational procedure to obtain a 
"real-time" identification algorithm that will run on the same time scale as 
the process to be controlled; (2) derivation of a convergence proof for the 
nonlinear least-squares identification algorithm that uses a "moving data 
window;" (3) generalization of the nonlinear dual-adaptive control method 
and derivation of asymptotic stability conditions for the resulting closed- 
loop system. 
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Abstract 

For larre complex structures , Knowledge needed 
for control system design cannot be determined ana- 
lytically at present. .Also, ground testing is not 
possible because of the sice of the objects and the 
fact that they are not designed to be self- 
supporting in Farth gravity. The aduptive/learning 
system for space operations assumes that structural 
testing must be conducted during deployment or 
assembly. Testing is conducted only when necessary 
to insure adequate control performance. Required 
desirn information is analytically extrapolated and 
monitored during other times. This report 
describes the adaptive/learning system proposed and 
presents simulation results using the SEP (solar 
electric propulsion) array and a novel remote 
sensor which involves raster scan television 
coverage of the motions of the array from four 
cameras on the corners of the space shuttle payload 
bay. This report includes a detailed description 
of the simulation, the filtering algorithm for 
processing the TV data, the parameter extraction 
algorithm, and simulation results gathered thusfar. 

I. Introduction 

Because of the desire for efficiency, the 
structural designs of large space systems are 
driven to low weight per unit area of surface. 

This necessarily results in low effective stiffness 
of the structures making them very flimsy and flex- 
ible. For such structures, operations involving 
barging, towing, and assembly require precise 
knowxedge of the structural dynamics of the objects 
involved to accomplish adequate attitude and con- 
figuration control. At present, for complex struc- 
tures, the knowledge required '-annot be determined 
analytically. Also, ground testing is not possible 
because of the size of the objects and the fact 
that they are not designed to be self-supporting in 
Earth gravity. Thus, it may be that an adaptive 
control system is required to deal with the prob- 
lem. Unfortunately, the control excitation levels 
usually required by adaptive systems (e.g. Ref. 1) 
are unacceptable for many precision space opera- 
tions. An alternate procedure is to schedule the 
adaptation process during those times when it is 
consistent with operational requirements and to 
analytically extrapolate and monitor the system 
during other times. An adaptive/learning system 
that operates in this manner was proposed and 
studied in an aircraft application (Ref. 2) and 
was later modified for use with structural systems 
(Ref. 3). The distinguishing feature of this 
system lies in the extrapolation and monitoring 
processes and the proposed system should not be 
confused with learning systems as described in 
reference 4 which deals with pattern classifica- 
tion. In the next section, the adaptive/learning 
concept is explained in the context of an example 
of the deployment of the solar electric propulsion 
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array (SEP) from the space shuttle- payload hay. 
Alternatives are discussed for the testing, •-xtra- 
■ Station, and monitoring subsystems. The estima- 
tion and parameter identification subsystems for 
which simulation results are avai’able are tner. 
mathematically outlined. This is followed ty i 
iescription of the simulation example and by 
rresentation and discussi' n simulation results. 


IT. A lai hive 7 learning ~;/.->tem Leper: : ti 

For the purpose of explaining the functions of 
the learning system, consider the problem of 
deploying a large solar array from its stowed con- 
figuration while attached to the space shuttle 
(Fig. 1). In the upper part of figure 2 a simpli- 
fied model is hown where the shuttle is repre- 
sented by the mass m and the solar array, fully 
deployed, in an undeflected state by the bar and a 
typical deflected state by the curved line. F 
represents the force resulting from the reaction 
Jets on the shuttle. The lower part of the figure 
shows two variations of the plots of the first 
fundamental frequency of the solar array as a 
function of the deployment length £ which is, in 
this case, the configuration variable. As indi- 
cated, the actual variation (shown dashed) might 
differ from that determined analytically. The 
system must initially use the analytic variation 
to scheuule the feedback gains during deployment. 

It is assumed that, at the start of the deployment, 
information is available through analysis or ground 
testing that adequately defines the dynamics, and, 
hence, produces adequate stability and control 
margins. In other words, the error between the 
actual and analytical frequency curves of figure 2 
is acceptable for small values of £. Also, if 
considered necessary, an in-flight test in the 
stowed configuration could be used to further 
reduce errors for 1=0. Semi-empirical methods 
must be used to incorporate the information 
obtained during the test into the analytic extra- 
polation used by the learning system. This is 
illustrated in figure 3 where the initial error 
(£ = 0) has been reduced (Test 1) and the analytic 
extrapolation has been adjusted as a result of 
Test 2. During further deployment, the error 
between the actual and the model frequencies may be 
expected to become large enough that further test- 
ing is required (test point 3). The decision to 
conduct another test can be put on an analytical 
basis by monitoring sensor signals during deploy- 
ment and comparing them to the output of an ana- 
lytic model based on the extrapolation. When the 
error grows to an unacceptable level, a new test is 
needed to insure favorable control system margins. 

For the adaptive/learning system discussed 
above, estimation and parameter identification 
theories are required together with an extrapola- 
tion of the system model, a method for evaluating 
the extrapolation to determine the need for test- 
ing, a method of designing the test inputs, and a 
method of refining the extrapolation process to 
bring it into agreement with the evolution of the 



.-.yjten. During tec tine : r trimeter identifiea- 
t ion a large signal t.o tioioe environment is 

rided - . . | teat inputs. In this -use , 

speed of testing is important because par-uneter 
identification of ultra low frequency systems - 
characteristic of largo space structures - will 
require considerable time. luring other ’ ines, no 
specific test inputs are .available. The only 
system inputs are those resuJ i.ig fron spacecraft 
perati.-ns. hence, there is i low signal to noise 
environment for parameter i ien. . fication. In this 
case, parameter identification in as low a sirnal 
to noise environment as possible is important 
assuming small deviations from the extrapolation 
model. For this case sequential parameter identi- 
fication as employed in reference 2 or a moving 
window .‘.'ewton-Paphson parameter identification may 
be desirable. The latter offers considerably more 
tolerance to noise but is also considerably more 
complex than the former. It also has monitoring 
information since estimates of the parameter 
variances are available by processing the informa- 
tion matrix of the data base (Refs. 5 and 6). The 
former is simpler to implement, requiring consid- 
erably les3 computation ani is amenable to parallel 
processing using a bank of second order parameter 
estimators. For that reason the former technique 
has been selected for development here. 

The technique used here was applied in 
reference 1 to the dynamics of a closed ring. The 
notion of the ring was assumed to be described by 
a Fourier series in the spatial variable with time 
varying coefficients. The coefficients were 
obtained using a fast-Fourier transform (FFT) of 
noise free measurements at sample times and their 
time series were used in parallel to identify 
frequency and damping characteristics of each 
coefficient. In that case a sine-cosine series was 
appropriate because of the closed nature of the 
ring. Generally, however, a sine-cosine approxima- 
tion may not yield a bank of uncoupled second 
order dynamic systems that describe the time evolu- 
tion of the coefficients of the FFT. In that case, 
the approximation functions of the series used to 
represent the dynamics must be "tuned" to effect 
decoupling of the bank of second order systems. 

The remainder of this section will deal with 
implementing the bank of sequential parameter 
estimators using arbitrary approximation functions 
and "tuning" the approximation functions to 
decouple the bank of parallel second order esti- 
mators. The general philosophy is to assume a set 
of approximation functions in the spatial and time 
domain (to account for rotating dynamics, e. g. 

Ref. 1), to obtain the "best fit" coefficients of 
a linear representation of the measurements at a 
sample time using the approximation functions, to 
select parameters of a bank of uncoupled linear 
second order difference equations to minimize the 
error between the model output and the observed 
coefficient time series, and finally, to "tune" the 
approximation functions to render the parameters 
time invariant. 

Analytically the motion, w, of the system can 
be represented by 


MM 

w(s,t,t) = l t|» (t,p U) )5 (s,c (0 ) (1) 

1=1 1 1 11 
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where C. ire the MM approximation functions to be 

used t > represent w; l i3 a measured configuration 
vs 

e . 


variable; are the time varying coefficients of 


.., is the snatial variable; t Is time; p. is the 

parameter vector of the i ' l linear second order 
difference equations modeling and = (c, , ) is 


, th 


ij 


the coefficient vector that relates the i func- 
tion £. to a constant set of M basis functions 

i 0. 


IJ 


(s,i) as follows: 
H 

’l 


I (s,c) = l a . ( s , i ) c . U) 
i j _ 0 IJ IJ 


(?) 


Thus, for each mode number, i, 
basis functions {o.,(s,0, J=0 


there is a set of 
* » »*^ ) which. 


when summed using the weights c. = (c. , j=0, 1,,, 

H ) according to equation (2), produces the i 
0 . 
i 

mode shape approximation function 5^. The c^ 


vectors will be adjusted, thus tuning the approxi- 
mation functions, only if monitoring indicates a 
deviation of the actual response from the extra- 
polation (l): One method for selecting the tuning 

functions o is illustrated in figure 1* where the 
* J 

mode shape is represented by a table of values 
at discrete values of s. The function is the 
best ground analysis estimate of the 1 st mode 


shape. Remaining a functions are selected to 
-LJ 

allow for adjustment of the value of 
corresponding to the j t ^ 1 table point. Initially, 


the values of +1 


dimensional c^ vectors are 


[l ,0, , ,0] which corresponds to using the best 
ground based analysis to start the system. 

Updating ^ and monitoring will require time series 


analysis of the parameters p... 


The parameters p.. are selected to minimize the 
error between the output of each model , 


^(k) = A li iJ/Jk-l) + A 2i ^(k-2) 

+ B n F i (k-1) + B a F i (k-2) (3) 


and the "best" fit of the sensor data to the 
approximation function representation ii.(k) at 

sample time t^. In (3) the parameter vector p. is 

p. = [A , A.., B, . , B-,, ] and F. is the time series 

force command inputs which are known to the 
controller. Thus iia is, in effect, the input 

sequence for a second order parameter estimator. 

In this report, this sequence is obtained using 
least squares fit of position measurements, y, 
which can be sequentially implemented. The form 
of the estimation problem is 
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NS 

y(k) = l H.(k)J.(k) ♦ v(k) (1*) 

i«l 


where y is the measurement vector, H. is a known 

operator algebraically related to C, v(k) is the 
error in the fit and J. is the estimate of ^ . 

The least squares solution to estimate ip is well 

known (Ref. 7) and will not be duplicated here. 

The solution nay be sequentially processed but one 
must take care not to propagate the estimated vari- 
ance of v. since this is equivalent to averaging 

the measurements at different times and will 
eventually produce an insensitive estimation of 1^. 

The technique used here to avoid this difficulty is 
to set the priori variance estimate to a constant 
at each k to start the sequential estimation. 


for all approximation functions, the identity e = 0 
(Eq. 5) hold3 for a constant p^ = p*. Consider 

then the variation of the IJM dimensional error 
vector F(p c ,c) = (e^) where e* is the i error 

function taken locally and p^ = (p c ^). The terms 

involving c appear explicitly in F only through the 
least square estimation process of i and hence we 
may calculate the change in F caused by a change in 
c by first calculating the change in ip given the 
measurement y (locally) and the change in c. This 
may be accomplished analytically (using the formu- 
las of Ref. 7) or numerically. Similarly a change 

in p results in a change in F which is linear 

c 

(through e) and is easily calculated. Hence, we 
may construct the equation 

iF = 7 F Ac ♦ V F ip 
c p r c 

c 


Subscripts i will now be dropped for simpli- 
city since processing for each approximation func- 
tion is parallel and identical. The parameter 
vector p can be selected iteratively using an out- 
put error formulation (Ref. 8) wherein the measure- 
ment sequence is (k indicating the sample time 

index) and the model output error is 

e k-l = V k-1 " Vk-2 ‘ V’k-3 


where 7 F = [V F,,,7F ] and Ap =[Ap 1 „Ap‘ ,M ] . 

p p p c c c 

c C 1 NM 

If a constant value p* = [ p . * ] exists then we may 

set Ap = p* - p , i = 1, 2 , , ,UM, and, hence, 
c i c ^ 

+ 7 F p = AF 

P„ c 


(V F, 7 F) 
C P c 


LP J 


- B l F k- 


2 " Vk-3 


(5) 


The iteration of p for the output error formulation 
is given by 


P k = p k-l + e k-l 


w E 

l v k-2 


W 2*k-3 


Vk-2 


W.F. o 
4 k-3 


( 6 ) 


In the last equation is the parameter vector of 

all modes identified locally assuming convergence 
of the output error algorithm. The solution of 
AF = 0 fcr p* and c provides, in the variational 
sense, the value of c that will result in a con- 
stant p* solution of the output error parameter 
identification problem for each mode taken over the 
data base. The next two sections concern the 
application of the methods presented here to a 
large space structure control problem. The only 
results gathered thusfar pertain to the least 
square estimation problem which serves as the input 
driver to the parallel output error identifiers and 
the output error identification process itself. 


As indicated in reference 1, the weights , ,W^ 
should be selected to be consistent with 


u i*L 


_2 


+ u I’ 1 " + W F‘" 

k-3 "3 k-2 Vk-3 


< 2 


and the inputs, 1 ^ and F , must be sufficiently 
varying and large if convergence of the sequency p 
to some value say p^ is to be assured. 

After convergence, the parameters track a 

sequence of converged parameters p which satisfy 

c 

(3) locally (for the sample sequence when the para- 
meters were identified.) If the sequence p^ is 

constant our .lob 1= done. However, if the sequence 
varies, "tuning" is needed for the approximation 
functions. The basic task is to adjust c so that, 


III. Simulation of the SEP Array 

A simulation of a large space structure has 
been developed to test the algorithms presented 
here as well as others. The structure used is the 
solar electric propulsion (SEP) array which is 
scheduled for deployment tests in orbit attached to 
the space shuttle as shown in figure 1. This 
section is a discussion of the finite element 
modeling of the array, the simulation, and the 
sensors and actuators simulated. 

The SEP array is illustrated in figure 5. It 
consists of a canister, an extendable mast, a con- 
tainment box with cover, and the solar array. The 
array is folded in the containment box and is 
lifted during deployment by the mast and box cover. 
The canister and box are attached to the shuttle 
payload bay using a carriage assembly not shown. 
First, the mast will be deployed without the array 
to test the mast deployment mechanism. Then, 
deployment of the complete configuration will be 
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tested. The deployment mechanism is designed so 
that it nay be stopped at of full extension or 
at full extension. At these points dynamics test- 
ing may be accomplished. He , we assume that such 
tests will be conducted by . *ring the space shuttle 
reaction control .lets which produce three axis 
force and moment control at the base of the 
canister. 

The simulation is developed as shown in 
figure 6. First, a finite element model of the 
physical structure is developed (Fig. 7) in a form 
required by the SPAR computer program (Ref. 9). 

The finite elerent model involves U 36 Joints with 
1705 interconnecting elements. The elements are 
bars, tubes, and beams which are included in the 
SPAR element repetoire A total of 1350 degrees of 
freedom ( including translation and rotational joint 
motion components as appropriate) are involved in 
the 3PAR model. A modal simulation was selected 
because of the large order of the state required 
to simulate the 1350 coupled degrees of freedom. 
Mode- frequency data required for the simulation is 
also provided by the SPAR program (Fig. 6). This 
mechanization also provides for the use of flight 
test data reduced to mode- frequency data (dashed 
blocks of Fig. 6). At present, a maximum of 10 
modes can be used in the simulation. However, that 
limitation can be removed without serious computa- 
tional penalty. 

The simulation also assumes a novel remote 
optical sensing concept that involves raster scan 
television coverage of the array (Fig. 8). Optical 
targets are distributed on the SEP structure and 
can be identified in the TV raster and registered 
in digital form at a given sample frequency. The 
sensors measure target motion but do not include 
components along the line of sight of each camera 
to the associated target. A total of 23 targets 
are assumed placed on the SEP structure which are 
observed with four cameras thus producing l81» 
sensor components at each sample time. The loca- 
tion of the cameras is shown in figure 8. (The 
targets on the containment box are not shown. ) All 
sensor data is generated during the simulation but, 
for the studies undertaken here, only l6 sensor 
camera-target combinations are used. These 16 have 
been selected to render some of the 10 modes 
essentially unobservable. 

In summary, the simulation to be used in 
evaluating algorithms for the adaptive/learning 
system is a modal one which uses up to 10 modes of 
six degree-of- freedom motion of the SEP array. The 
simulation inputs are the forces and moments 
applied to the base of the SEP canister. Sensor 
information is the motion of spots located on the 
SEP structure as perceived by four comeras located 
in the space shuttle payload bay. 


IV. Simulation Results 

The simulation results gathered thusfar per- 
tain to tne first test — mast deployment without the 
attached array — at 3 A and full extension. The 
least squares estimation of modal amplitudes has 
been developed and tested, in simulation, assuming 
that o.q functions are the ones computed from SPAR. 
The sensor targets were three equally spaced along 
the containment box cover (at the top of the mast) 
and one spot halfway down the extension axis on the 
mast. All four cameras were used. This spot- 


camera combination renders the seventh SPAR mode 
essentially unobservable since it is prelominuntly 
a torsional excitation of the containment box, the 
targets on which were not processed. The assumed 
a ? r ior i variance of 5 was proportional to the 
square of the product of frequency of the associ- 
ated mode and the sample time, t = .1 sec. For the 
array simulated the lowest model frequency was 
.109 Hz and the highest frequency (mode 10) was 

Hz. (This will obviously challenge the sample 
data parameter identification process). The modes 
are ordered in increasing frequency and Shannon's 
sampling theorem is violated for modes 8, 9, and 
10. Mode 7 is, as mentioned earlier, unobservable. 
Also, modes 1 and 2 are very close, frequency-wise, 
but correspond to vibrations which are totally 
different spatially. Measurement noise was 
assumed to be independent of other measurements 
and Gaussian with zero mean. Studies have been 
made with measurement noise standard deviation of 
.03 in. to .003 in. dependent on the range of the 
target to the camera. For the results, herein, the 
sensor standard deviation was .01 in. for the 
closest target which was on the mast and .003 in. 
for the targets on the containment box cover. 

Thus, the sensor processing was essentially noise 
free with respect to motions of about six inches. 

Figure 9 shows the spatial distribution of the 
mast motion along its deployment axis at selected 
times during a 10 sec. simulation at 3/ 1 * extension. 
Only the x^ component is shown which is normal to 

the space shuttle longitudinal axis for 6=0 
(Fig. 8). The simulation resulted from a unit 
excitation of all 10 modes. 

The measurements taken from cameras 1 and 2 
(x^ and X 2 raster components) for the target on the 
mast are shown in figure 10. Figures 11, 12, and 
13 are results taken from the estimation and 
identification algorithms of the learning system 
during the simulation. Figure 11 corresponds to 
controller mode 2, figure 12 to mode 6, and 
figure 13 to mode 10. In each figure we have the 
sequence of 5! estimates and the identifiers itera- 
tive estimate of A^ and A^. Figure 13 also 

includes the modal estimation error which was the 
largest of all 10 modes. The length of time for 
convergence is substantially larger for mode 2 
parameters than for mode 6. This is, of course, to 
be expected since the frequency of mode 2 is r 
lower than that of 6. Mode 7 results are not 
presented since cannot be estimated. For this 

mode the parameter identification algorithm does 
not converge using the sequence as input since 

it is uncorrelated with the motion. Tests as 
recommended in reference 5 should be made on the 
estimated variance of ij before passing it to the 

sequential identification algorithm. The mode 10 
results are somewhat unexpected since convergence 
to the actual parameters does occur although 
Shannon's sampling theorem is violated. Table 1 
shows the actual, initial, and final parameter 

values for all modes. The B parameter values were 
not identified since there was no force or moment 
inputs. Note that under the assumed structure the 
unobservability of mode 7 does not affect other 
modes. 
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This taper has presented an adaptive/learning 
Control system concept for large space structures. 
The concept employs :t priori ground tests and 
analysis initially to r.odel the structure. The 
learning process requires in-flight testing to 
refine the model whenever control system perfor- 
mance degrades to an unacceptable level. Hub- 
systems required for learning involve distributed 
sensor processing and parameter identification — 
algorithms which have been proposed in this report. 
A simulation of a large space structure has been 
developed to study the performance of the proposed 
algorit.ams. The structure simulated was the solar 
electric propulsion array which is scheduled to be 
flight tested for deployment and structural 
dynamics. To this time the only identification 
performed has been the parameters of a bank of 
second order difference equations (one for each 
mode in the learning system model). For the 
simulations conducted the parameter estimator 
behaved as expected — converging to the actual 
parameter values. The behavior of the system with 
realistic noise levels in the sensors is yet to be 
evaluated as is the approximation functions 
"tuning" algorithm which was presented. 
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Table 1 Initial, Final, and Actual A^ 
Parameter Values for the Simulation 


Mode 

Initial 

Final 

Actual 

1 

.9953 

1.9953 

1.9953 

2 

.9952 

1.9902 

1.9952 

3 

.9855 

1.9852 

1.9855 

L 

.31*70 

1.31*1*3 

1.31*70 

5 

.2637 

1.2637 

1.2637 

6 

-1.51*82 

-.51*82 

-.51*82 

7* 

-1.981*9 

I .8561 

-.981*9 

8 

-2.81*23 

-1.81*12 

-1.81*23 

9 

-2.9973 

-2.0291 

-1.9973 

10 

.7520 

1.71*13 

1.7520 


‘Unobservable 
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Figure 1.- Solar Array Deployment from 
the Space Shuttle. 




I 

Figure 2. -'Simplified model of solar array deploy- 
ment showing the variation of frequency 
with array length. 
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Figure 8.- 'HP Remote Sensing Concept. 
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Figure 9.- Mast Deflection Along the 
Deployment Axis. 
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Abstrac 


Algorithms are presented for on-line parameter 
identification cf structural dynamic systems. As 
an example, they are osed to calculate the para- 
meters of a modal model of a flexible beam. The 
algorithms are tested using hardware consisting of 
a 12 ft. beam with four voice coil actuators and 
r.ir.e ncr.ccntacting displacement sensors. They are 
cror rmre d ir. a CSC Cyber l”5i iigital computer 
which provides input command signals for the act- 
uators, reads the sensor data, and processes .the 
algorithm to calculate consistent estimates of the 
modal parameters of the beam. Experimental results 
are compared with those of simulation analysis. 




r. 


The control of larw 
will depend cn knowledge 
of the objects involved, 
this information can be 
from ground testing, ait 
been proposed £--2] • 
ir indirect ittro’ic'n tc 
ir. the clock diagram of 
based upon a modal decor. 
spor.se of the flexible : 


e flexible space structures 
of the structural dynamics 
Since it is unlikely that 
determined analytically or 
poive control schemes have 
Tr.e of these schemes [ 3] , 
adaptive control, is shewn 
Figure 1. This approach is 
ipcsition of the dynamic re- 
tructure ar.d was designed 


to make use *f the parallel processing features of 
modem microcomputers. Tr.e of the significant 
problems involved in the implementation of this in- 
direct approach to tr.e adaptive control of large 
space structures is that of real-time, or on-lir.e, 
identification of the parameters of the large, flex- 
ible structural components. Satisfactory perfor- 
mance of the parallel structure identification tech- 
nique indicated in Figure 1 car. be achieved only 
when the appreximatien functions noted in the Figure 
correspond to the natural modes of the flex.b'.e 
structure £o] . This paper presents a new st unique 
for using cn-lir.e measurements to improve estimates 
of the approximation functions, ind combines this 
approximation function up- dating apprcach with the 
medal parameter identification algorithm creser.ted 
in [3]'. 


Time series analysis for the purpose of para- 
meter identification for multi-input, multi-output 
systems has received a reel leal of attention C' — ”] 
in recent years under various assumptions retarding 
the parameter matrices appearin' in these models 
ar.d cn the nature of assumed random disturbances. 


reference [•*] contains a thorough treatment of 
least-squares and instrumental viriabie approaches 
to the problem of generating consistent estimates 
of parameters in univariate an i multi variate models 
The instrumental- vari a ole approach to system ident- 
ification is examined in [ o j for a single-input , 

* Professor, electrical Engineering 

**Aero-3pace Technologist, Member AIAA. 


single-output system where -easurer.er.ts cf 


sys- 


tem input and noisy measurements of the system out- 
put are ava.lable. It is shown following [-] that 
tr.e instrumental variable method yields consistent 
caroneter estimates. In £ o ^ a review of a number 
of methods of estimation of parameters ir. multi- 
variate autoregressive moving average models is 
creser.ted. Among the techniques ccr.s i dered are 


wer.eralit 
methods , 
cedures. 
proposed 
shown to 
identical 


ed least-squares, instrumental variable 
and maximum likelihood estimation crc- 
A "limited information" estimation is 
in [6] and, for autoregressive models, is 
produce parameter estimates asymptotically 
to conditional maximum likelihood estimat- 


es while requiring reduce i computational effort. 

In C a statistically efficient least-squares pro- 
cedure is leveled e d for estimating the natural 
frequencies and damping parameters cf a system 
characterized by a lumped mass-spring- camper model 
under stationary random forcing functions. The 
only available measurement is a sequence of iispl- 
aoements subject to additive r.oise. Although the 
formulation in does not include the treatment 
of known forcing functions r.cr the inclusion of 
multi-output measurements, numerical results pre- 
sented demonstrate that the derived least-squares 
estimates are a good approximation to maximum-like- 
lihood estimates. 


In the following section the structural para- 
meter identification problem is formulates as a 
problem of calculating the parameters of a linear 
discrete-time multi-input , multi-output system 
using r.oisy on-lir.e measurements. Ir. section 3 we 
combine the identification algorithm of Co] with 
the new approach for updating the approximation 
functions. The MASA-Langley structural dynamics 
and control test facility is then described ir. 
Section 1. Experimental results presented in 
Section 5 are obtainei by applying the pr:posei on- 
line identification proced'ure to a flexible beam 
assuming that the approximation functions are those 
taken from an a priori simulation or from a ground 
test of the flexible structure. A comparison with 
the results of digital computer simulation studies 
is also included. 

— . Prcbl 4 ?!*'. f c r’~ u 1 1 * i c n 


The sccpe cf the identification problem car. be 
examined fren the following formulation: the 

motion v(s,t) of a flexible structure may be re- 
cresented by 

v* 


vher* s ier.ctes Vr.e 
ser.ts time. 3-i - s * 
tior.s, for example. 


spatial ••ariable and t r 
finite set of If. spatial 
the -'irst :!M eigenfur.cti 



,s ir. 


a modal expansion, are time-varying coeffi- 
cients of , or moial amplitudes, ar.d v ( s , ti 

retresents these higher-order terms r.ot appearing 
in the finite sum. The parameter vector 


1 


represents a set of parameters ir. a difference 
equation governing the time-evolution of the 
fi ( t h > Px) , :'3r example 




ani 

. For later 
as 


where * ( ** ^ T, ~k > Pi ) 

*>i = CA,;,A lXl B,;,B,J 

reierer.oe we oesigr.ate the ccmccr.er.ts of fc- 

: ■ V"i 3 : 'PJ ~ A 'i ' F* = Azi , pjs 

K =i? > i • l- tr.ere are '.In actuators applying 
t.-.e set of farce inputs U i , A = / > . . . > :IA, 
to the flexible structures, then the r.oial forces 
Fg at tearing in v 2) are giver, by 


VA 


“F(o =y >7,.. U(ii) 

j * / ^ ; 


(3) 


The measurements available to the 
controller are a set of sample i values of the 
notion at the points , j = 1,2,.. ./J S- 

NM 

^ (l,)= Il (S j)?. (1 «.fc) tV( W*^<!) (U1 

4= / 

where HS is the number of ser.scr locations ar.i 


’J 


•eser.ts the measurement errer 


Sr 

n 

U,(w" 

?w= 

’v* 

F0*)= 

fib) 

U(kk 

Mlk) 

• 

• 

• 


7.1 s ' 




F(k) 

** J 




If we write the set of measurements as the 
vector ^(k) , the set of modal amplitudes as the 
vector |(nj, the set of actuator ani modal forces 
as the -ectors U(k) and F(k) , respectively. 


(5) 


* v* have t r. e fo llo vln«k node! 

= Nj(k) + W(k) (6) 

The identification problem i3 the on-line calcul- 
aticr. of the elements of the matrices appearing 
in >;), where A (< A Z , S, , B, ire iiagcr.aj mat- 
ri:es with tair. diagonal ei.ener.ts A (i , A, i( 3 M( 3^; 
t = /, , . . , IIM respectively, " is an X-! x ilA 
matrix and H i3 an ii3 /. til! matrix whose columns 
are linearly independent in the model model. The 
error vector v(k) appearing in (6) is an U’S-vector 
whose components le/",- (fc) are the sum of the measure- 
ment errors e t (V) ar.d the moiel errors vtS tj r^) 

U' e i (fcU'/'(S t ,t # ) j x- l t . . .NS ,7! 

If we define , then the model !o) 

simplifies to 
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*(*> = -*■ u^fe) 


(3) 


The available measurements are the series of act- 
uator forces f ar.d sensor measurements 

{*(*>} • 


"wT.— pop* p *■. ® f 1 d? 


- i f 1 c it i c n 


The adaptive control system. in figure 1 was 
designed to tahe advantage of the parallel pro- 
cessing capability cf m.oder microcomputers, 'dr. for- 
tunately, the parallel architecture is correct only 
in oase the approximation functions ir. the figure 
are the natural nodes of the struct -rsl system, 
cecause of tr.e imser.se computational problem posed 
by computing the moce shapes, it is unrealistic 
to assume that the control outputs, pcle-tero 
tarameters, ar.d me :e shapes car. be computes at the 
same . requer.cy. the philouptiy adopted herein is to 
assume that the control outputs are computed at the 
highest frequency, parameters that define the pole- 
zero characteristics at a lower frequency, and the 
cose shapes are updated when possible. This sec- 
tion first summarises ar. algorithm, ta/.en from 
reference^], for identifying the parameters that 
define the pole-zero characteristics of the system. 
Then an original method for updating the approx- 
imation functions in Figure 1 is presented. 


1. Identification of the ?:le-Cero Character'. st i os 

Here it assumed that the approximation functions 
are given and that the parallel architecture of 
figure 1 is valid. Under these assumptions the 
problem is reduced to one of identifying 

of equation (2). The 

technique of reference I ol is followed here. It 
is ah output error formulation wherein the error 
is given by 




(9 1 


ani the formula for updating the estimates of the 
vector n? is 


£(10 = k(fc-i) + ev(A-i) 


ir.g ' 

the 

! es 

t i’. 

\ 

b 

(k- 

2) 


ii 

(fe- 

3) 

W i 

K 

(k- 

2) 

w. 

I* 

(h* 

3) 


(13) 


In equation (1C) the weights « should be selected 
:o satisfy the requirement that 


W, £ ( b- z)* Wj q*(k- i) * VJ } f M f W [Ik-]) <2 5 

to insure stability of the identification algorithm. 

2. U’ndatir.w of the Ancrcximation Functions 

The updating approach is based upon rewriting 
the second equation in ,3) as 



2i*)= H *(h) 
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£y] r rt T {kr"i) 


H*= (fc'wjV 


mo . • ier.otes ratrix transpose. The inverse 
ippearir.? in \ 1 2 ' exists sir.ce the colu.rr .3 of X 
ire linearly inieper.oer.t. rest lain? ( 12 ) with 
the first equation ir. (3) yielis 




r N 

ij(k> n*i)= m/ f n(l;*N*l) 

n, t 

»I. 

Equation (IT) nay he written ir. r.ore tort act 


Y{u) = S(N)M T + V(n) 


M< =MA,M , *./,* 


, < = 1,1 ( 15 b) 

-n (te ♦ l) = »r(k'l)- M, »'(k)-M^\*-(k- l) CWe) 


S(N)=[Y(N-)iY(M- 2 )|u(N-.)jir(N- 2 )] 


< = 1,1 

( 15 o ) 

Y\k) = 

j T (»<-(N-3)) 

0 

uo= 

U(L-(N-3)j 






u \ l * 2 ) 

• « 


K = k ♦ M - bsN-2 L = R + N- 


It is interest in? tc note that if )J3 a )!M, ie. 
ve have selected a r.ur.cer of measurements equal 
tr.e trier of the r.oial rcdel to he iier.tifiei, 
ir. )15a) r.ay he -.Titter, as 


M it»* *[>> « 


r^ r (fe+i) 



y Vx/ # 


n r (fe*Wt 2 ) 

O 

o 

G 

r 


, is a iia- 


iener.ts^ij-r 

"cte that the Y(K) 

A is ar. ei?er.- 

(J(L) are ilxi.'A mat 


M 

. K 

M " N. r 


Her.oe, l preoeiure for ieterr.inir.g the nodal 
pararecero is as follows: 1) Use neasurer.er.t iata 

to iete.--.ine the nitrites M, , ‘I ^ i n . 11 ). 

3 liltulate the eigenvalues and eigenvectors of 
:!,ar.i i ' 2 to iete—ir.e A,, A., ana W . 2) halve 

fir 3 ( mi 3 2 fron .lot). Since the eigenvectors 
of i* A ire linearly independent hut are, in .-er.eni, 
non-unique, the iier.tifitatior. of H mu hence, gf 
2, mi u 2 is i.so non-unique. This, however, is 
hot surprisir.? riven the fom of the nodal model 


The tatrioes ir. :) .-.ay he ieterr.ir.ee ty tah- 
ir.i i sequence of -.exsureser.t iata *.(*“). y 
• . .jj^OfWtJ) lhi orierir.e these measurements 
is follows: 


s'oo = 


* 

e 

^ (Kt 2 ) 


Ucw, if + NA) , a least-squares solu- 

tion M T to (if) nay he obtained as follows: 

M T = S Y(n) ( 22 ) 

v-ere S^ier.ot es a tseuio-inverse that is readily 

f:,puted [5.3j 


3 
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However, lines the disturbances V(>V,)are con- 
temporaneously correlated with the eler.er.ts of 
5(N)the ordinary least-squares procedure will 
yield asym ptotically biased estimates [•] . A 

technique for providing consistent estimates is 
inj - run^nt 2l viriibl^ T.^thoi ^ ^ . Z f 

the la,-«e i output and input measurements j 

V,W.O(k-0,-.., ^^•0.%(^2\V< k,N - , ),V(k*N-s) 

are causer. as tar. uii-.te instrumental variaties, 
we fern the array £ r ( /v ) , 


Z(N) = 


123 ; 



• • j ( k *■ n~ i) 


. . }(Jf» A/- 2) 

\J(k) • 

. . U(krN-l) 

V(k-) • 

. . If(ktM-Z) 


rron the iefiniticn of the instrumental variable 
matrix Z (M) in '22. ar.i measurement error matrix 
V in (20) and (15) it is tlear that the cor- 
responding columns of Zf*7ani rows of V(N) are 
unoorrelated. 

r 

If we premultiply (It) byZ(fj)ve obtain 

Z (MY(n ) -2 (n)$(n)N) -t-Z(N)V(N) ' 21, 

” - a ir.s t r’^jr.snt al r.is * r.*? i*om 

M r = [zmSm] £ r MY{ n) ( 2 5) 


A computational procedure based upon the 
foregoing discussion is as follows: 

1. Assume the model to be identified is of order 
‘2! ar.i, from a set of measurements , form the 
N*2(H AeNM) matrix SlN)in (le). 

2. Form the instrumental variable matrix 7(w)ind 
the solution (25!* 


3. 

ar.v) 


Compute the eigenvalues of the ; X 

the eigenvectors of the /w|. , v- v*. ( = ) 2 

Zcmrute the solution B ) , to 


A A , 

hB J = n j 


4 . ) 


i =1,2 


(26) 


7 rcr. these steps the i i e r. 1 1 !* i e i rtcial rr.o del 
for to ie j 1.3 


The solution 25) can also be expressed in an 
iterative fashion [5,61. 

It is interesting to note that tr.il procedure 
also provides an independent estimate of the modal 
pole-oero parameters. However, because the f cre- 
ating procedure requires the computation of a 
matrix pseudo-inverse ar.d eigenvectors ar.d eigen- 
v a lues of the computed matrices M, and 1'^ it is 
anticipated that the up rating of the approximation 
functions will take place at a lower frequency 
than the computation cf the modal parameters from 
the output error algorithm. 


IV. Experimental Test Eacilitv 


The algorithms oreser.ted above h :. e beer, 
tested by digital computer simulation and by ap- 
plication to the UASA-Lar.gley structural dynamics 
ar.i control experimental apparatus. 

The experimental test facility 

is shown in the photograph of Figure 2 and in the 
diagram of Figure 3. The apparatus comprises a 
uniform rectagular-section 12 ft. long flexible 
beam, four electromagnetic actuators, nine r.cn- 
ccntacting sensors to measure deflection cf the 
team at various locations, ar.d four strain gage 
load ceils at the actuator attachment points. 
.Although the sensors ar.i actuators car. be placed 
at arbitrary stations along the beam, the locations 
of the sensors ar.i actuators used in the exper- 
iments reported below are those shewn sor.e.-.atical- 
ly ir. Figure 5. 

The beam is supported using two cables that 
are attached to the building. These are used to 
off load the bearings of the four actuators showr. 
in the photograph of Figure 2. The actuators are 
’Jnhoit c-Diokie Vcuel 1 Shakers which have a stroke 
of 1" and a maximum force output of JO Ibf. They 
are attached to the beam with spring steel flexure 
couplers that permit transverse vibrations in the 
horizontal plane while suppressing torsional mo- 
tions. The couples are attached to strain gage 
load cells which were manufactured in-house at 
La?C. The load cells are then hard mounted to the 
team. Ar. additional purpose of the cables is to 
maintain a horizontal loading condition by pre- 
venting iroop of the end of the actuator attach- 
ment assembly. Cn the side of the beam opposite 
the actuators are noncontacting deflection sensors. 
These are Kam.ar. KT-23CC series system sensors. 

The range of the Hamor. sensors is 2 in. The 
resolution of the sensors is *.CC1 ir. . 

The components mentioned above are interfaced 
to the laPC Tyter 175 real time computer system, 
as indicated in Figure The Cyber computer 
reads the sensor data, processes the algorithm to 
calculate estimates of the parameters of the beam, 
ar.d provides input command signals for the act- 
uators. A high gair. analog electrical circuit is 
used to assure that command forces from tr.e con- 
trol system, resident or. the Cyber computer are 
realized at the lead cells on the team. For that 
reason the force commands from the Cyber computer 
are sent to the EAI 650 analog computer where the 
oircuit adjusts the current delivered to the act- 
uators so that the commanded load occurs at the 
load sells. The maximum csrmar.dei load from the 
Cyber is set tc 1C Ibf whereas the maximum value 


ORIGINAL PAGE fs 

OF POOR QUALITY 




5 



O *jT« 



ORIGINAL PAGE !S 
OF POOR QUALITY 


Figure 2.- Experimental beam facility. 
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modes of the beam. 
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Figure 12.- Output of parameter Identifiers, mode 3 excited (experimental) 
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Igure 13.- Output of pararaet 
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Abstract 


This paper presents the formulation of sim- 
ultaneous state and parameter estimation problems 
for flexible structures in terms of least-squares 
minimization problems. The approach combines an 
on-line order determination algorithm, with 
least-squares algorithms for finding estimates of 
modal approximation functions, modal amplitudes, 
and modal parameters. The approach combines pre- 
vious results on separable nonlinear least squares 
estimation with a regression analysis formulation 
of the state estimation problem. The technique 
makes use of sequential Householder transforma- 
tions. This allows for sequential accumulation 
of matrices required during the identification 
process. The technique is used to identify the 
modal parameters of a flexible beam. 

1 . 1 nt roduc t i on 

On line identification of parameters and 
states of distributed parameter systems is a 
significant problem affecting the design of 
adaptive controllers that have been proposed for 
large flexible space structures. in [ l] we 
proposed a technique for identifying the modal 
approximation functions and modal parameters of a 
f 1 ni t e-d imons i ona 1 approximate model of a flexi- 
ble structure. The technique proposed in [ lj 
was based upon a priori knowledge of NM, the num- 
ber ol modes being excited during an identifica- 
tion time interval. Because the number of excited 
modes fs generally unknown we propose here a 
procedure for identifying NM on-line. This tech- 
nique is based upon using Householder transforma- 
tions [2] to determine the rank of a sequential- 
ly obtained measurement matrix. Once NM is iden- 
tified, the approximation function updating pro- 
cedure of [ ij Is used to determine a set of 
normalized approximation functions. The modal 
parameters are then determined using a non-linear 
least-squares parameter and state estimation al- 
gorithm. This algorithm, reported in [3j , is 
based upon the work of Kaufman [ij and Golub and 
Pcreyra [ 3] on separable nonlinear least squares 
prob 1 ems . 

In the lol lowing section the structural 
parameter Identification problem is formulated as 
a problem o( estimating the parameters of a linear 
multi-input, mu 1 1 i -on tput system of unknown order 
using nofsv on-line measurements. The identifi- 
cation algorithms are presented in section 3 and 


applied in section U to identify the parameters 
of a flexible beam from measurements performed at 
the NASA-Langley structural dynamics and control 
test facility. 


II. Problem Formulation 


The motion w(s,t) of a flexible structure 
is represented by 

w-(S,-t) = O, (t fc )jjt) + Vfi.ti (1> 
i=' 0 

where s denotes the spatial variable and t repre- 
sents time. The represent the first 

NM eigenfunctions in a modal expansion, 

) are ">odal amplitudes, and 

v(s,U represents those higher-order terms not 
appearing in the finite sum. The parameter vec- 
tor represents a set of parameters charac- 

terizing each mode. For example, li each mode is 
characterized by the differential equation 

'L. + 2: u q. .+ zoffii) ( 2 ) 

[A J. I t’A *0* t J A' 

then parameter vector would comprise the com- 
ponents ( J",’ / , A',‘ J . If each mode were 

characterized by a difference equation 


f. ■(**')= A, ,1 (k)rA^(k-i)* Bffft) r jF(bi) (3) 


then parameter vector % f\ ^ would comprise the 

components ( A,;, A Jw s ( -, 3 • in ( 2 ) 

and (3) -j x ft ? and 3=*.^' respectively represent 
the modal forces in the continuous and discrete 
mode 1 s . 


The measurements available to the controller 
are a set of sampled values of the motion at the 
points s.^ , j = 1 , 2 , . . . ,NS , and a set of sampled 

values of the actuator forces \J A , 1*1,2 NA, 

where we have assumed NS sensors and NA actuators 
are used. Thus, if Jf (/c! denotes a measurement of 
the flexible structure motion at point s- , then 


1 {*) =1^1 )j 4 i f >, h , ) + mV*) t?(k) < 4 > 


where 



represents the measurement error and 
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force applied 
t , then 


If Ktk) represents the modal 
to the i-lh flexible mode at time 

HA 

T U}» n- [7 (fc) 

j" 


(5) 


where i=l,2,...,NM and where the r> Aj actuator in- 
fluence parameters are usually unknown. 


The parameter identification problem to be 
treated below is that of on-line processing of the 
available sequences of actuator force measurements 
{Ul(fc),} and sensors measurements 

in order to estimate the number of excited modes 
NM, the approximation functions the 

modal parameter vectors , and the actuator 

influence parameters • . 

111. Identification Algorithms 


A set of identification algorithms Is des- 
cribed below for performing on-line parameter i- 
dent i f ication in the following sequence: (1) est- 
imate , . M .,(2)estimate the J (S,).and( 3 >est ima t e the 
^andmjV From the description below it will be 
seen that while the computational procedures for 
accomplishing these tasks make use of common vec- 
tors of stored measurement data, the computation- 
al time for each procedure differs from the others 
so that, in general, these identification subtasks 
will operate in parallel with transfer of updated 
estimates occuring as new estimates become avail- 
able . 

1. Identification of number of exc ■ led modes NM 


To examine the nature of this entification 
subproblem, define the vector of measurements U(kj 

- 


and modal amplitudes 


by 



' *(U* 

. ?w* 


* 

* 

H ik) 


S? (k) 

. 


Ohm 

. 


( 6 ) 


so that from (A) we have 


^ 6 ?) = H 


(7) 


where the NSxNM matrix H has jij -th element 



and the measurement error vector w(k) is an NS- 
vector whose components are given by 


i'k) * V(K,t k ) 


(9) 


If we denote the columns of 11 by h - * ( 

J*1 ,2, . . . ,NM, then c7) may be rewritten as 


^ (k) - 4 + ...+ £ jjk) t ur(k) ( 10 * 


In (10) NM as well as the vectors h>) are un- 
known and must be determined from the on-line 
measurements y(k). The procedure oescribed below 
for finding NM is a generalization and extension 
of the work of Woodside to multi-dimension- 

al systems. 


First consider the case where w(k) = 0. The 
problem then is to determine the number of lin- 
early independent vectors hi required to repre- 
sent the measured response vectors y(k)=z(k), 
where 


NM 

}(*)= X 

j*' 

Assume that there is some maximum system order 
NM" < NS beyond which we will not attempt to 
model the flexible structure. Form the array 


( 11 ) 



From (11) it is seen that ifM is such that 
NM" > M > NM, the M columns of jJ-(M)will be 
linearly dependent whereas, the columns of-j2(ttf 
be linearly independent if M £ NM. Hence, to 
determine NM we examine the rank of the^-(M) 
array for a sequence of M values. If it is found 
that there is a value M* such that 

ran ra * feljlMVOj (13) 


then NM=M*. 


A simple rank test is based on the orthogonal 
decomposition of the data matrix (12). As shown 
in [2j if A is an mxn matrix (m > n) of rank 
k<>\ then there is an mxm orthogonal matrix 
Q and an nxn permutation matrix P such that 
it > 


Q A P = 



<lt 

O 1 


t 


m. 


ij 


(1A) 


where "R ,, is a fexk upper triangular matrix of 
rank h . A computational algorithm based upon 
use of Householder - rans f orma t ions for obtaining 
the decomposition (1*») is specified in [2] and 
is applied in section A below where the data 


3 


matrix ^(M) replaces matrix A. Testing the main 
diagonal elements ol the right hand side of (14) 
provides a rank test and, hence, an estimate of 
the number of excited modes. 

2 , Identification of Modal Approximation Func- 
tions 


Once the number of excited modes NM has been 
determined, a least-squares procedure can be 
used to estimate the modal approximation func- 
tions. A procedure was presented in £ l] for de- 
termining the modal approximation functions 
■V)... a discrete modal description of the 

form 


^(kti) =A,*lb)-* A.yft-O* B.'J(k) ♦B.U'Jr-i) ( 1 s ) 


with available measurements (7). In (15) A ( and 
A^ are diagonal matrices, q(fc ) is a vector of 
modal amplitudes and IK J? ) is a vector of act- 
uator forces. The procedure consists in finding 
the eigenvectors corresponding to the NM largest 
eigenvalues of estimates of matrices M ( and Mj 
where 

if (fe ♦<)= M, ^ (0 f M, ij (*-()♦■ NU(k)+fy 1 b ) 


and 

KrHAH* } K (17) 


yi (n a i} = l)-M, * »r(k)~ Mj ‘* r A’- /) ( 18 ) 


H * in (17) is given by 

H* = (H t H)''H t < 19 > 


where 


r r 1 I | 1 

S = I V (r) -i^ • Y (*-2^ > >)■ \J 'N -2 1 |(21) 

l i i . - 


"rM= 

/(K- 

• 

1 

U(L)= 

^ {l-(h-i)) ! 

• 



Of 

f- 

J 

V 

K= fc »n ( Imw-i , 

R.N-2 

L- 

J 

and 

v'ffiij) 


M. T " 

\ 

i 

v(n) = 

• 

« 

• 

m t - 

Mi 

N. r 



f T ln*^t2) 





(22) 


(23) 


From (17) it is seen that the columns of H are 
eigenvectors of M^ and the non-zero eigenvalues 
of M^ are the diagonal elements of the matrices 
A ^ . Furthermore the matrices can he deter- 

mined from (17) once estimates of H and are 
available. This procedure requires the storage 
of measurement data over a time-interval that in- 
cludes at least 2(NM-NA) measurement samples. In 
situations where more rapid up-dating ol parameter 
estimates is desired, alternate techniques £ll , 

[ 3J may be used. These are described below. 

3. Identification of Modal Parameters 

In [l"! it was assumed that the estimates of 
modal approximation functions would be calculated 
at a lower frequency than the computation of the 
modal parameter estimates. The latter were ob- 
tained from an output error algorithm. In this 
paper we employ a nonlinear least-squares algor- 
ithm [3J which could also be used to provide 
estimates of modal position and velocity. The 
technique is based upon computing a least -squares 
estimate of modal amplitudes at each sampling 
time. For each mode i we use a state description 
of the form 


where ( • ) denotes matrix transpose. 

Least-squares estimates of M,, , N ( , and 

Njare obtained by finding the least-squares 
matrix solution M to 


»*(k*i)r F'* 4 (k)+ Gr 4 (fc)L T (0 < 24 > 

where 


( 20 ) 


A 

h 


O I 


Mn)= SCn)M t + V (N ) 


(25) 
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and C* is a 2x5 matrix whose elements depend on 
the modal parameters vector ^ . From this 
point on we will drop the superscripts and sub- 
scripts "i" since the following algorithmic app- 
roach is to be used for each of the NM excited 
modes . 

Since we assume here that the modal approxi- 
mation functions comprising the H matrix of (7) 
have already been identified, we now use (7) at 
each sampling instant to generate an estimate 

q< L ). 


^{O = *, (k) ♦ r>(k) 

- [/,o] x(k) f -»(*) 


(26) 


of modal pos i t ion ,whe re n ( k ) denotes the error 
In the modal amplitude estimate. A separable 
nonlinear lea st -squa res formulation of the par- 
ameter/state estimation problem is obtained as 
follows: Solve (24) with respect to U ( k ) to 

obtain 


U(k) - 6 f (r>*M- (/(rJRrJ^fO 


-GhiRr' G%;] 

’ r( k) ' 

U IP ll'; M { > J 



where C (p) denotes the pseudo-inverse of G(p). 
Combine (26) and (27) as suggested by Duncan and 
Horn [ 7"j to form the overdetermined system 


u (ol 
3 o i 
V(D\ 

NU)'= F(r) 

' ^ e- 

=. Z • ' 

t 

r 
+ u- 

r n 

A. . 

'n f !) 

yni) 

• 

h • 

# 

K 

• 

Uf Kj 



»(*) 

• J 


i(*j 

r j 

where 




I • 

H 

■ (tirWf), 0 Vr) ••• 
H • •• 


and 


C- = 


. . . Off) 

... H 


(29) 


( c 
I C 1 


( 30 ) 


The vector on the left-hand side ol ( 26 ) is a set 
of input-output measurements tawen over an in- 
terval of length K. Lquation (26) may be rewrit- 
ten more briefly as 


Y(k> F( r )X(K) + G„ e(K) 


(31 ) 


The simultaneous state/paramet cr estimation pro- 
blem is now considered to be the least-squares 
minimization, 


|| V(K)jj‘ (32) 


Following the work of [dj and ' 5 ' , (32) is 
separated into the form 


Vni n { || RXM-QtfiOlj 2 * j Q, Y(h')f f (33) 


where Q, , Q 2 , and R are defined from the orth- 
ogonal decomposition of F k (p), 



R ! 

o\ 


(34) 


It is shown in [_ 3 , 5 J that the optimum parameter 
estimate p and optimum state trajectory X(K ) 
that mininize (33) can be obtained as solutions 
to two subproblems: a nonlinear least-squares 
problem 


>r»in || Q , (f) V(K)||^ ( 35 ) 

to yield p , and a linear 1 ea st-squa re s problem 
based upon 

ft(p) X(k,£)s- Q,(p) (36) 


to yield X(K). 

An algorithm for the solution of these 
problems based upon the approach of Kaufman 
has been implemented and used in the simulation 
study reported in the following section. 
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IV. Simulation Results 

The algorithms described above have been 
implemented m a simulation program for identi- 
fying the characteristics of a flexible beam 
used in the NASA-Langley structural dynamics 
and control test facility. The experimental test 
facility is described in detail in £l] . Six 
flexible modes were used in simulating the beam 
response. Figure 1 shows the results of apply- 
ing the algorithm of section 3 for testing the 
number of excited modes when only three modes 
were actually excited in the simulation. In the 
figure the magnitude of the main diagonal element 

of the test matrix on the right hand side of 
(14) is plotted against the number of assumed 
excited modes. Fesults of tests with noiseless 
measurement data and with two levels of noisy 
data are presented. It is seen that when the 
assumed number of modes exceeds the actual num- 
ber of excited modes <33 the corresponding mam 
diagonal element remains below the test level 

T 

Figure 2 shows the result of applying the 
regression analysis approach of section 3 to 
identifying the mode shape functions when it was 
established that only one flexible mode was 
excited. Test results for two cases are shown 
in the Figure: (a) a simulation wherein the beam 
was given an initial displacement corresponding 
to the first flexible mode and (b) a laboratory 
experiment in which only the first flexible mode 
was excited. In the simulation the normalized 
mode shape obtained agreed with that produced by 
the SPAR structural analvsis program. In the 
laboratory experiment a distortion in the ident- 
fied mode shape was observed. In both cases the 
mode shape was identified on-line after approx- 
imately one second of data processing. 

Figure 3 presents one illustration to dem- 
onstrate the convergence of the nonlinear least 
squares parameter identification approach of 
section 3 for identifying the modal parameters of 
the first flexible mode. Four samples of meas- 
urement data are used and convergence of the 
algorithm is achieved after ten Iterations. A 
study of the algorithm performarct with various 
levels of measurement noise is presented in ^3^ . 

V . Cone lus l on 

In this paper the st a t e/pa rametcr identifi- 
cation problem for large flexible systems is dec- 
omposed into a number of simpler subproblems 
which may be solved by least-squares algorithms. 
The technique comprises the following steps: 

(1) estimating the number of excited modes in a 
set of sequential measurement data, (2) estima- 
ting the spatial approximation functions needed 
In a modal response representation and A 3 ) estim- 
ating the pole-zero parameters and moon 1 position 
and velocity coordinates lor each excited mode. 

Particular least-squares algorithms for ac- 
complishing the three subtasks above were imple- 
mented in a simulation for Identifying the para- 
meters of a flexible beam. These algorithms 
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employ Householder transformations thus avoiding 
the need for matrix inversion and resulting in 
accurate and rapid identification. Other techni- 
ques mav also be considered tor solving each sub- 
problem. For example, in testing the number of 
excited modes, identification of the me* sui ement 
noise covariance matrix could lead to a more ac- 
curate setting of the test level used. One sub- 
ject of current research is a comparison of com- 
peting techniques with regard to computational 
requirements and achievable accuracy. 
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Figure 1. lest for Number of Excited Modes 
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ABSTRACT 

A technique is presented for obtaining a control law to 
regulate the nodal dynamics and identify the modal parameters 
of a flexible structure. The method is based on using a 
min-max performance index to derive a control law which may 
be considered to be a best compromise between optimum 
one-step control and identification inputs. Features of the 
approach are demonstrated by a computer simulation of the 
controlled modal response of a flexible beam. 

I. INTRODUCTION 

A class of indirect adaptive control systems proposed for the 
control of large space structures [1] is based on a modal 
decomposition of the system dynamics and may incorporate one 
or more on-line testing schemes [2] to determine when 
successful parameter identification has been achieved. The 
control strategy used in calculating the actuator inputs must 
achieve adequate regulation or tracking performance and, at 
the same time, provide inputs to allow adequate parameter 
identification. A control system designer is thus faced with 
Llie problem of devising a control strategy to ensure 
acceptable system performance even when on-line parameter 
i denti f i abi 1 i ty tests have failed because the system 
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configuration has changed or the environment in which the 
system operates has changed. 

In this paper we formulate and examine the performance of a 
nonlinear dual -adapti ve control scheme in which a 
sampled-data controller is designed to select a best 
compromise between an input signal that is optimum for 
mean-square system regulation and an input signal that is 
optimum for parameter identification. Dual control theory, 
originally formulated by Feldbaum [3,4], has been studied in 
[5-7] and in the references cited therein. A key concept 
introduced by Feldbaum is the dual control strategy based on 
a performance index that takes into account the fact that 
future observations on the process will be made. A 
controller may be able to "probe" the system for state and 
parameter estimation improvement, which then may improve 
future regulation and tracking performance. In many 
situations where the dual nature of stochastic control is no: 
taken into account the controller becomes "cautious" [5,6] 
and tends to "turn-off". This undesirable phenomenon is 
avoided by the approach described below. 

II. FORMULATION OF AN ADAPTIVE PERFORMANCE INDEX 


The di sc rete- ti me dynamics for each mode is assumed to be 
described by the ARM A model 


y ( t ) +aiy ( t- 1 ) + a2y ( t-2 ) = bi u ( t- 1 ) + b 2 U ( t-2 ) +e ( t ) (1) 


where y(t) denotes modal displacement, u ( t ) denotes modal 
force, and e(t) is a sequence of independent, 
equally-distributed, normal ( 0 , o 2 ) random variables. It is 
assumed that e(t) is independent of y ( t- 1 ) ,y ( t-2 ) , . . . , 
u( t-1) ,u( t-2) , . . . and that the parameters a 1 .a 2 .b 1 .b 2 
are unknown constants. If we let Yt denote the information 
available to the controller at time t, 

Yt = jy ( t ) ,y ( t- 1 ) u( t-1) ,u( t-2) , . . . j. . (2) 

x(t) denote the modal parameter vector and e(t) denote a 
modal measurement vector, 


x T (t) = (ai ,a 2 f bi,b 2 ) ; 

e T (t) = ( -y( t-1 ) ,-y( t-2) ,u( t-1) ,u( t-2) 

where ( . ) T denotes vector or matrix transpose, then (1) 
may be rewritten as 

y ( t ) = o T ( t ) x ( t)+e( t) 


(3) 


(4) 



where the constant parameter "dynamics" satisfies 


x ( t + 1 ) = x ( t ) (5) 

It can than be shown, following the analysis of [8], that the 
conditional distribution of x(t + 2) given Y t+ i is normal 
with mean x(t + 2) and covariance matrix P(t + 2) where x(t) and 
P(t) satisfies the difference equations 

x(fl) = x(t)*K(t)(y(t)-eT(t)x(t)) (6) 

K ( t ) = P(t)e(t)/(o2+ e T(t)P(t)©(t)) (7) 

P(t+1) = P(t)-(P(t)e{ t)eT( t)P(t) )/ 

( o2+®T ( t ) P ( t ) © ( t ) ) (8) 


Furthermore, the control law that minimizes the regulation 
c r i t e r i o n 

V c ( u ( t ) ) = E-[y 2 (tn) 1 Y t | 
is given by 

y‘Vi.(t*l)x,(t*l)+P, 4 (t*l))e,(t*l) 

u(t) = - i 

x| ( t"*- 1 ) + P 33 ( t+1) 

where ^ denotes the sum over i = 1 to 4 with the value 3 
exc 1 uded. 


(9) 


( 10 ) 


To provide bounded modal inputs that improve parameter 
identification accuracy while guaranteeing that the modal 
amplitude will not become excessively large, the controller is 
designed to optimize, at each sampling instant t, the following 
performance criterion: 

min max [ V ( a , u ( t ) ) ] (11) 

u( t) a 

subject to the constraints 

um^M.O^A^l (12) 

where 


V c (u(t)) 


V j ( u ( t ) ) 


v ( a , uun 


A 


c 


(1-X) 


(13) 



V c denotes an acceptable or desired level of regulation cost. 
Vj(u(t)) denotes and identification cost function of u(t), 

V I( u( t) ) = trace ^P(t + 2)J (14) 

e 

Vj denotes and acceptable or desired level of identification 

cost. The maximization indicated in (11) yields a function 
V ( u ( t ) ) which, although not convex, is interpreted as 
specifying, for each admissible u(t), the most costly linear 
combination of relative regulation and relative identification 
cost. Minimization of V ( u ) thus yields the modal input that 
minimizes this most costly combination of relative 
identification and regulation performance. 

III. SIMULATION RESULTS 

Since V c ( u ( t ) ) and trace P(t + 2) are relatively simple 
functions of u(t) the numerical solution of the one-step 
optimization problem (11)-(13) at each sampling time is quite 
feasible. Results of simulation studies described below 
illustrate an interesting feature of this approach: since the 

parameters involved in the evaluation of V c (u(t)) and 
V j ( u ( t ) ) depend on system measu r|ments , the optimum 
distribution of relative cost, x°(u) depends on on-line 
measurement data and hence, at each sampling instant, the 
weighting between identification and regulation will change 
depending on the on-line system performance. This is in 
contrast to[9] in which a fixed weighting between absolute 
control and identification cost is used at each sample time. 

In the simulation study we compare the performance of three 
control systems: 

a) A constrained adaptive controller that minimizes (9) 
subject to the control magnitude constraint. 

b) An optimum identification controller that minimizes 
(14) subject to the control magnitude constraint. 

c) The one-step dual-adaptive controller based on 
( 1 1 ) - ( 13 ) . 

In Figures 1-3 we present simulated modal response data for the 
first flexible mode of the Langley beam experiment described in 
[10] where we assume here that a single actuator is used. The 
accumulated on-line regulation cost, VT, shown in Figure 1 is 
defined as 

■ & 


VT ( N ) 


(15) 



and the on-line identification cost, PT, is defined as 


PT ( N ) = trace [P(N)] (16) 

where P ( N ) is calculated on-line using (8). Note that for the 
first 10 to 15 sampling times the regulation cost of the 
dual - adapti ve controller is close to that of the constrained 
mi nimum-vari ance controller and the identification cost of the 
dual - adapti ve control system is close to that of the 
constrained one-step optimum identification controller. 

Figure 2 indicates that the dual -adapti ve controller's actuator 
signals switch between its limits, _M).5, more frequently than 
do the actuator signals of the other controllers. This may be 
due to the lack of any energy constraint in the above problem 
f ormul ati on . 

A future study will examine the performance of the 
energy-constrai ned dual -adapti ve controller in comparison witn 
energy-constra i ned mi nimum-vari ance and one-step optimum 
identification controllers. The relative regulation cost and 
relative identification cost defined in (13) are plotted in 
Figure 3 where 

V c (N) = o 2 N (17) 

is the accumulated control cost that would be achieved if the 
parameters of the system where known precisely and if an 
unconstrained control law were used; o2 = 10-4 was used in 
the simulation runs. A constant value Vj r = 10-4 was chosen 
as indicating the acceptable level of parameter 
identification. Figure 3 indicates that, depending on on-line 
measurements, the one-step identification and regulation cost 
at one sampling instant can have widely differing shapes from 
their respective distributions at other sampling times. This 
leads to the on-line variations in the dual -adapti ve control 
strategy mentioned earlier. 

The simulation results indication that the one-step, 
constrained dual -adapti ve controller has the feature of 
providing, based on measured data, system inputs that result in 
parameter identification while maintaining bounded modal 
amplitude response. 
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